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LONG-TERM  GOALS 

The  goal  of  this  research  is  to  acquire  a  quantitative  understanding,  leading  to  predictive  models,  of  the 
broader  aspects  of  linear  and  nonlinear  sound  scattering,  transmission  and  coherency  in  oceanic  bubbly 
mixtures  pertinent  to  the  shallow-water  acoustics.  This  includes  a  conceptual  understanding  of  the  role 
played  by  damping  and  stabilization  mechanisms  in  bubble  dynamics  and  longevity.  Of  particular 
interest  is  the  delineation  of  different  regimes  of  behavior  as  a  function  of  frequency,  size  distribution, 
flow  and  volume  fraction. 

OBJECTIVES 

An  objective  specific  to  this  project  is  the  extension  of  the  theory  of  sound  transmission  and  coherency 
in  bubbly  liquids  to  derive  attenuation  characteristics  for  both  small  amplitude  (linear  response)  and 
large  amplitude  (nonlinear  response)  forcing,  ultimately  incorporating  the  effects  of  contaminating 
surface-active  solutes.  A  second  objective  is  the  development  of  a  unique  laboratory  capability  for  the 
precise  and  accurate  measurement  of  the  frequency-dependent  complex  acoustic  impedance  of, 
scattering  from,  and  the  coherency  of  propagation  through  well-characterized  bubble  ensembles  for 
frequencies  spanning  the  individual  bubble  resonance  frequencies.  Characterization  implies  the 
precise  knowledge  of  the  space-  and  time-dependent  bubble  density  and  size  distribution. 

APPROACH 

The  approach  involves  a  balance  between  theory,  analytical  modeling  and  experiments  to  predict  and 
measure  propagation,  scattering,  and  coherency  characteristics.  The  dynamics  of  a  single  bubble  for 
both  small  and  large  amplitude  forcing  and  is  treated  numerically  using  the  Keller  formulation  for 
bubble  dynamics.  Attributes  of  bubble  behavior  (mainly  damping  and  resonance  response)  can  be 
quantified  and  incorporated  into  a  comprehensive  description  of  sound  propagation,  scattering  and 
coherency  by  extending  the  Wood-Foldy -Morse  theories.  The  final  step  is  to  incorporate  the  effect  of 
surface  active  materials  by  adapting  numerical  models  developed  by  Church  [1995]  and  Allen  &  Roy 
[2000a&b],  among  others. 
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Laboratory  experiments  eover  two  fronts  of  aetivity.  FYOO  efforts  foeused  on  the  measurement  of  the 
eomplex  impedanee  of  bubble  distributions  terminating  a  sound-hard  impedanee  tube  over  frequeneies 
ranging  from  well  below  to  well  above  bubble  resonanee.  The  bubbly  medium  ean  be  eharaeterized 
optieally  using  a  stereo  mieroseope  and  eleotrieally  by  way  of  eonduetivity  measurements  to  determine 
void  fraetion.  A  variant  of  the  dual-sensor  impedanee  tube  teehnique  is  employed.  From  this,  the 
frequeney-dependent  phase  veloeity  and  attenuation  of  the  bubbly  medium  is  obtained. 

The  seeond  experimental  thrust  for  FYOO  is  the  eonduet  of  tank  experiments  to  measure  baekseatter 
from,  and  transmission  through,  bubble-filled  test  seetions.  Both  spherieal  and  eylindrieal  geometries 
are  employed.  Experiments  are  eondueted  at  frequeneies  higher  than  4  kHz  and  are  to  supplement  thin 
walled  tube  experiments  at  frequeneies  between  15  and  30  kHz.  The  targets  employ  bubbly  media  for 
whieh  the  eharaeteristies  are  known  and  boundary  eonditions  that  are  rigorously  defined.  These 
experiments  are  also  a  preeursor  to  the  FY02  field  experiments. 

Modeling  is  eurrently  foeused  on  the  development  of  a  eomprehensive  and  self  eontained  formalism 
for  deseribing  the  aeousties  properties  of  bubbly  liquids,  drawing  from  the  elassieal  efforts  of  Wood, 
Foldy,  Morse  and  others.  This  is  supplemented  by  a  parallel  effort  aimed  at  developing  Ist-order 
eorreetion  terms  that  aeeount  for  the  surfaee  elastieity  and  dilatational  viseosity  of  surfaee-aetive  layers 
similar  to  those  that  ean  eoat  oeeanie  bubbles. 

WORK  COMPLETED 

A  eentral  thrust  of  these  aetivities  is  to  eomplete  the  development  of  experimental  teehniques  suitable 
for  use  in  field  test  planned  during  the  FY02.  During  FYOO,  an  ONR  "Foeused  Aeousties  Bubble 
Workshop"  was  held  at  Boston  Univ.  to  determine  the  researeh  issues  of  oeeanie  bubbles  in  shallow 
water  and  their  influenee  on  propagation.  On  behalf  of  the  ONR,  we  invited  experts  in  this  field  to 
determine  the  seope  of  these  issues  and  to  develop  a  researeh  plan  that  provides  for  their  resolution. 
Workshop  partieipants  diseussed  the  major  theoretieal  and  experimental  issues  as  well  as  the  type  of 
shallow  water  experiments  to  be  eondueted  in  eoordination  with  the  High  Frequeney  Aeousties 
Communieation  program  during  FY02. 


In  FYOO,  preliminary  laboratory  experiments  were  eondueted  to  explore  sound  seattering  from  well- 
defined  assemblages  of  bubbly  fluid.  To  date,  two  primary  advanees  have  been  made.  First,  a  bubble 
generation  system  has  been  developed  to  produee  large  amounts  of  appropriately  sized  bubbly  fluid, 
from  whieh  to  make  the  experimental  targets.  Seeond,  an  initial  round  of  seattering  experiments  has 
been  eondueted  using  targets  made  with  this  new  bubble  generator.  Targets  also  ineluded  voided 
polyurethane  spheres  and  thin-walled  polyurethane  &  polyethylene  spheres  filled  with  bubbly  gel. 

At  low  frequeneies,  bubble  elouds  ean  be  modeled,  to  first  order,  as  aeoustieally  eompaet  eompressible 
spheres.  An  effeetive  medium  approximation  [Wood,  1930]  is  used  to  infer  the  bubble  eloud's 
eompressibility  based  on  the  void  fraetion.  This  was  demonstrated  at  Lake  Seneea  [Roy  et  al.,  1992] 
using  artifieially  generated,  freely  rising  bubble  elouds.  Present  experiments  are  being  eondueted  in  an 
indoor  tank  using  smaller  bubbles  with  higher  resonanee  frequeneies  in  order  to  span  the  effeetive 
medium,  resonanee,  and  the  sonie  veloeimeter  regions.  Several  liters  of  bubbly  fluid  are  neeessary  to 
eonstruet  the  desired  seattering  targets  and,  until  now,  there  has  been  no  way  to  generate  these  types  of 
fluids.  Shown  in  Fig.  1  is  a  sehematie  of  a  new  high-volume  bubble-generation  deviee.  Water  is 


pumped  through  a  porous  ceramic  pipe  is  impregnated  by  compressed  air  forced  through  the  wall  of 
the  pipe,  resulting  in  a  nearly  mono-disperse  distribution  of  bubble  sizes.  The  sizes  are  inversely 
proportional  to  the  liquid  velocity;  hence  a  degree  of  control  is  available.  Recirculation  is  utilized  to 
increase  the  void  fraction.  Samples  are  tapped  from  the  main  processing  system,  pass  through  the 
scattering  test  section  (not  shown)  and  returned  to  the  holding  tank.  The  flow  rate  through  this 
secondary  circuit  is  low  enough  to  be  considered  static  on  the  time  scale  of  a  single  acoustic  pulse. 


Figure  1:  (a)  Schematic  of  the  bubble  generation  system.  (b)Typical  bubble  size  distribution 
produced  by  the  system.  The  sample  was  taken  from  a  test  section  at  STP. 


Initial  scattering  experiments  were  conducted  at  the  NUWC  Acoustic  Test  Facility  using  voided 
polyurethane  spheres,  bubbly-gel  fdled  spheres,  and  a  10-m  test  section  of  thin  walled  latex  tubing 
through  which  bubbly  fluid  flows.  For  the  initial  analysis,  the  targets  were  modeled  as  uniform,  wall¬ 
less  compressible  spheres  and  cylinders  (infinite  length,  see  [Stanton,  1988]).  Independent  void 
fraction  determination  is  not  yet  available,  so  VF  was  fit  as  the  only  adjustable  model  parameter. 
Scattering  results  from  cylinders  with  three  different  void  fractions  are  shown  in  Fig.  2.  A  good 
quantitative  fit  is  obtained  between  theory  and  measurement  for  the  two  lowest  void  fractions,  but 
there  is  a  6  dB  difference  for  the  highest  VF  case.  Although  our  initial  analysis  does  not  account  for 
this  discrepancy,  there  nonetheless  am  intriguing  qualitative  match.  The  resolution  of  this  discrepancy 
calls  for  further  experimentation/analysis  however,  the  initial  results  are  most  encouraging. 


Figure  2:  Measured  and  predicted  values  of  scattering  from  a  bubbly-fluid-filled  latex  tube  with 
1.27-cm  id.  Solid  lines  represent  model  predictions  using  void  fraction  to  fit  the  measurements. 
Below  8  kHz,  the  scatter  in  the  two  lowest  VF  cases  may  be  caused  by  poor  signal-to-noise  ratio. 

The  impedance  tube  apparatus  was  developed  (see  Fig.  3)  to  characterize  the  bubbly  liquids  and  initial 
measurements  were  reported  last  year.  Results  agree  with  model  calculations  in  all  respects  except  for 
frequency-dependent  damping.  During  FYOO  we  have  formulated  a  damping  model  with  accurately 
accounts  for  radiation  damping  and  will  be  incorporated  into  the  propagation  models  for  both  the 
thick-walled  impedance  tube  and  the  thin- walled  pulse  tube.  Impedance  tube  results  were  also 
obtained  using  a  bubble-free  polymer  termination  and  a  water  termination  (i.e.  an  open  ended  tube).  In 
both  cases  the  data  agreed  with  theoretical  predictions,  however  artifacts  were  observed  in  the  data  that 
suggest  that  a  more  accurate  sensor  phase  calibration  is  required  when  probing  samples  with  low 
intrinsic  attenuation.  An  improved  calibration  scheme  is  nearing  completion. 
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Figure  3.  Impedance  tube  apparatus  used  to  characterize  bubbly  fluids. 


There  were  two  other  techniques  employed  for  generating  bubbles  for  laboratory  scattering  and 
propagation  experiments.  Freely  rising  micro-bubbles  were  generated  using  a  rotating  porous  cylinder 
though  which  air  is  forced  [Chiba  &  Takahashi,  1998].  Bubbly  distributions  were  also  created  by 
suspending  micro-bubbles  in  a  1%  solution  of  Xanthan  gum.  The  advantage  of  the  latter  technique  is 
the  bubble  distributions  in  both  size  and  space  are  stationary  in  time  and  easy  to  characterize.  The 
disadvantage  is  that  the  impact  of  the  viscoelastic  suspending  gel  is  not  known.  To  address  this  latter 
point,  we  intend  to  measure  the  acoustic  response  of  a  single  bubble  suspended  in  the  polymer  gel. 

The  frequency-dependent  response  of  the  bubble  will  be  measured  using  laser  Mie  scattering  and 
compared  with  model  predictions  of  bubble  dynamics  in  Newtonian  and  viscoelastic  media. 

Since  the  recognition  that  oceanic  microbubbles  have  a  dramatic  influence  on  sound  production, 
scattering  and  transmission,  many  articles  have  appeared  in  the  acoustics  literature  concerning  the 
properties  of  these  bubbly  mixtures.  To  extend  this  work  to  turbulent  media  we  have  examined  the 
engineering  literature  on  two-phase  fluid  flow.  Relevant  experimental  results  have  been  summarized 
and  a  unified  theoretical  treatment  has  been  prepared  to  provide  clarity  to  the  community  as  well  as  a 
basis  for  our  experimental  work.  This  includes  an  extension  to  the  Morse  and  Ingard  treatment  for  the 
scattering  and  transmission  of  sound  through  turbulence  containing  bubbles.  This  treatment  provides 
the  basis  for  the  calculation  of  signal  coherency  and  it  is  our  intent  to  test  this  theory  in  careful  tank 
experiments.  Moreover,  we  have  made  considerable  progress  in  developing  a  model  of  sound 
propagation  through  bubbles  coated  with  contaminating  skins  and  layers.  Correction  terms  that 
account  for  surface  elasticity  and  dilatational  viscosity  have  been  incorporated  into  the  force  balance 
equation  [Jankovsky  et  al,  1999]. 

RESULTS 

Multi- frequency  backscattering  experiments  have  been  performed  in  a  tank  using  a  variety  of  bubbly 
targets  of  known  geometries.  Voided  polyurethane  sphere  target  strengths  (2-20kHz)  were  measured 
for  void  fractions  of  0%,  3.4%,  and  6%;  target  strengths  for  the  voided  spheres  where  on  the  order  of 
40  dB  {re  Im).  The  frequency  response  exhibited  modal  structure,  with  peaks  shifting  to  lower 
frequencies  for  lower  void  fractions.  No  backscatter  signal  was  detected  for  the  solid  polyurethane 
sphere.  Target  strength  was  also  measured  for  a  hollow  polyurethane  and  polyethylene  spheres 
containing  a  suspension  of  bubbles  in  polymer  gel,  and  from  a  bubbly-fiuid-filled  latex  tub.  Results 
are  consistent  with  an  effective  fluid  scattering  theory.  These  initial  results  are  preliminary  as 
independent  void  fraction  determination  measurements  were  not  performed.  However,  qualitative 
agreement  was  found  with  reasonable  estimates.  The  latex  tube  experiment  used  our  new  method  to 
generate  large  volumes  of  mono-dispersed  bubbly  fluid  samples. 

Model  calculations  for  surfactant  coated  bubbles  show  an  increase  in  attenuation  for  frequencies  far 
removed  from  the  bubble  resonance  for  radii  less  than  100  microns.  There  is  no  significant  increase  in 
attenuation  for  larger  bubbles  (order  1  mm).  Interestingly,  near  resonance,  the  interfacial  viscous  case 
has  a  lower  attenuation  than  the  clean  interface,  most  likely  due  to  the  reduced  radiation  damping 
associated  with  smaller  radial  excursions.  Likewise,  resonance  effects  on  the  phase  speed  dispersion 
curve  are  also  significantly  diminished  due  to  the  interfacial  viscosity. 


IMPACT/APPLICATIONS 


The  notion  that  bubbles  ean  be  driven  to  pulsate  eolleetive  is  important  to  any  assessment  of  seattering 
and  attenuation  from  oeeanie  bubble  elouds  and  layers.  The  area  of  researeh  is  important  to  HF/SW 
noise  and  propagation,  SW  mine  hunting  sonars,  high  power  aeoustie  arrays  for  MCM,  and  wake 
homing  torpedoes.  Furthermore,  the  aeoustieal  measurement  of  bubble  populations  and  eireulation 
patterns  may  depend  on  the  physies  of  multiple  seattering  and  absorption  in  bubbly  mixtures.  Finally, 
the  persistenee  of  mierobubbles  in  the  shallow  water  eolumn  ean  also  effeet  the  resolution  of  sonars 
and  these  effeets  need  to  be  quantified. 

RELATED  PROJECTS 

I .  A  eollaboration  with  NUWC/Newport  (R.  Costa)  and  CSS/Panama  City  (K.  Commander)  on  an 
instrument  designed  for  in  situ  measurement  of  oeeanie  bubble  dissolution  rates.  The  instrument 
was  sueeessfully  deployed  as  part  of  a  wake  physies  sea  test  at  Nanoose  in  September,  1998  and 
the  results  were  published  in  the  Proeeedings  of  Oeeans  ’99  [Costa  et  ah,  1999]. 

2  .  A  eollaboration  with  APL/UW  (S.  Kargl)  on  the  physies  of  nonlinear  beam  forming.  The  issues  of 

bubble-enhaneed  nonlinearity  and  dissipation  in  bubbly  water  are  key  to  both  projeets.  The  goal  of 
the  beam  forming  projeet  is  to  assess  the  viability  of  high-intensity  sound  beams  for  MCM. 

3  .  A  eollaboration  with  Univ.  Virginia  (J.  Allen)  on  developing  new  theories  that  deseribe  the 

dynamies  of  bubbles  in  viseoelastie  media.  The  PI  was  the  researeh  advisor  for  J.  Allen,  a  UW 
Ph.D.  student  in  Meehanieal  Engineering  who  eompleted  his  dissertation  September  ‘97. 

4.  A  eollabration  with  the  NRL  has  been  proposed.  Meetings  have  been  held  with  Dr.  Franehii  and  his 
staff  at  NRL. 
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